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ABSTRACT 

This report documents the full details of the condensed journal article by Ashpis k 
Reshotko (JFM, 1990) entitled ’’The Vibrating Ribbon Problem Revisited”. A revised 
formal solution of the vibrating ribbon problem of hydrodynamic stability is presented. 
The initial formulation of Gaster (JFM, 1965) is modified by application of the Briggs 
method and a careful treatment of the complex double Fourier transform inversions. 
Expressions are obtained in a natural way for the discrete spectrum as well as for the 
four branches of the continuous spectra. These correspond to discrete and branch-cut 
singularities in the complex wave-number plane. The solutions from the continuous 
spectra decay both upstream and downstream of the ribbon, with the decay in the 
upstream direction being much more rapid than that in the downstream direction. 
Comments and clarification of related prior work are made. 
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Preface 


This report documents the full details of the concise article entitled "The Vibrating 
Ribbon Problem Revisted” by Ashpis & Reshotko which appeared in the Journal of Fluid 
Mechanics (1990), vol. 213, pp. 531-547. This report constitutes the full manuscript that 
was originally submitted for publication, but had to be shortened and condensed due to 
space limitations. It is based on Ashpis & Reshotko (1986), a work principally supported 
by the US Air Force Office of Scientific Research. The present report was written after the 
first author joined NASA. Some typographical errors of the prior work were corrected and 
comments of the journal reviewers were implemented. There has been recently renewed 
interest in the continuous spectrum in context of its relationship to the new concepts of 
transient growth in hydrodynamic stability theory and its role in receptivity to freestream 
disturbances. Therefore it was felt that the details and the structure of this report make 
its documentation in form of a NASA report a valuable complement to the journal article 
and an accessible text to non-specialists. 
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1. Introduction 


The vibrating ribbon became a common experimental device in hydrodynamic stability 
research since its first use by Schubauer and Skramstad (1947) to excite TS waves in 
boundary layer flow. The importance of this fundamental device goes beyond this 
practical application. The vibrating ribbon problem serves as a simple example of 
receptivity problems, which express the effect of imposed disturbances on the flow, and 
are significant in describing the onset of instabilities leading to transition from laminar 
to turbulent flow. 

This problem, also known as the signaling problem, was initially treated by Gaster 
(1965), whose analysis verified Schubauer and Skramstad’s observation that in the time- 
asymptotic (long time) limit the ribbon excites the spatial eigenmodes of the flow at its 
own frequency. Additional analytical work is reported in Russian literature (e.g. Tumin 
and Federov, 1984). In free shear layers the signaling problem was treated by Huerre k 
Monkewitz (1985) for the inviscid case. We found Gaster’s solution incomplete. Mainly 
missing are the continuous spectra. Therefore we provide here a corrected and complete 
solution. 

Gaster’s initial formulation of the problem is followed: the model, the formulation 
as an IB VP (Initial Boundary Value Problem), and the application of a double Fourier 
transform (or Laplace-Fourier transform), are essentially similar, but we proceed with 
a different solution procedure. We use the Briggs method (Briggs, 1964) to obtain the 
time-asymptotic solution. Common for long time in plasma physics (see review by Bers, 
1983), this method was first used in fluid mechanics by Tam (1971, 1978). It became 
more established in this field following a large number of its applications in later years. 
For example it was used by Huerre k Monkewitz (1985), Ashpis k Reshotko (1985), 
Leib k Goldstein (1986), Pierrehumbert (1986), Monkewitz k Sohn (1986), Hultgren 
k Aggarwal (1987), Monkewitz (1988), Lin k Lian (1989), Yang k Zebib (1989), and 
was reviewed by Huerre (1987). The Briggs method requires an accurate account of the 
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singularities in the transformed planes which is performed here in a methodical way. As 
a result the continuous spectra are included in our solution, and various related topics 
are illuminated. 

The continuous spectra were addressed by Case (1960, 1961) t. A compehensive 
treatment was given by Grosch Sc Salwen (1978) and Salwen Sc Grosch (1981), who 
extended the classical normal mode formulation. They identified the one branch of 
continuous spectra in the temporal case, and the four branches in the spatial case, but left 
unsettled a difficulty in the physical interpretation of two of the branches in the spatial 
case. 

Continuous spectra were obtained also in formulations of linear stability problems as 
I VP (Initial Value Problem) and IB VP: When solved by Fourier or Laplace Transform 
methods, the continuous spectra emerge from branch-cuts in the complex transform 
planes. The temporal continuous spectrum was obtained by Gustavsson (1979) who 
formulated an IVP. His work was used by Salwen Sc Grosch (1981) to prove completeness 
of the temporal eigenfunction expansion. Tsuge Sc Rogler (1983) treated a forced problem 
formulated as IBVP and obtained two branches of the spatial continuous spectra, but 
dismissed one of them as nonphysical. They ignored an inconsistency in their solution, 
which yielded a branch-cut singularity in a region of the complex transformed plane that 
should be analytic. The same branches were obtained also by Aldoss (1982). A similar 
forced problem for compressible flow was also treated as an IBVP by Tumin Sc Fedorov 
(1983), who found seven branches of the spatial continuous spectra. They recognize the 
same inconsistency, proposed a method to resolve the difficulty, but which does not seem 
to be a practical one. 

In the solution here presented, the continuous spectra emerge in a natural way as 
part of the complete solution and the solution is mathematically consistent. The solution 

t After completion of this manuscript the authors have found that the continuous spectra was also 
addressed by Murdock &: Stewartson (1977) and Hultgren & Gustavsson (1981). 
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procedure is applicable to a wider class of problems, and clarifies the prior difficulties. 

The outline of this work is as follows: The problem is formulated in Section 2 and 
the solution in the transformed domain is presented in Section 3. Reference is made to 
necessary details which are given in Appendices A and B. The inversion to the physical 
domain is performed in Section 4, which includes description of the singularities and 
the application of the Briggs method. The physical meaning of the results is described 
in Section 5. The discussion in Section 6 includes a comparison to Gaster (1965), and 
comments about related topics. It is found necessary to include an outline of Gaster’s 
solution in Appendix C in order to simplify the comparison. The present work is based on 
unpublished work by the authors, Ashpis & Reshotko (1985,1986), where further detail 

can be found. 
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2. Formulation of the problem 


The physical problem is modeled as shown in Figure 1. Two-dimensional 
incompressible boundary layer flow over a flat plate is assumed, and the ribbon is 
modeled as a line source of disturbances imbedded at the wall at the origin O of the 
cartesian coordinates x, y. The longitudinal and normal velocity components are U and 
V, respectively, and the freestream velocity is Uoo • We are interested in the response 
of the boundary layer to two-dimensional harmoni'c excitation of the ribbon starting at 
t = 0. 

The assumption is that the vibrations of the ribbon are infinitesimally small, 
justifying use of linear stability theory. The velocity is split into basic ( U,V ) and 
disturbed (u, v) parts, and the standard procedures of linear stability theory are applied. 
Assuming parallel flow, V = 0, U = U(y) , the nondimensional equation for the normal 
disturbance velocity v is 


8 Tr ® 1 r -?2 ~ tt " d v 1 *-,4 

m +u T* Vv ~ u T* = H Vv 


( 1 ) 


where R is the Reynolds number based on Uoo and on the displacement thickness 8* , 

{?' = — , V 2 = — — -f —— , and v = v(y:x,t). The independent variables are 
dy 2 ox 2 dy 2 

used in the order y; x,t throughout this paper for convenience. 

The effect of the ribbon is formulated as a wall boundary condition on the normal velocity 


u(0;x,t) = cos wo t 8(x) H(t ) 


(2a) 


where u>q is the real frequency of vibration of the ribbon, 8 is Dirac’s Delta function, 
which expresses the modelling of the ribbon as a line disturbance, and H(t ) is the unit 
step function, which indicates that the motion starts from rest at t = 0. The second 
boundary condition at the wall, from continuity, is 


dr(0;x,t) 

dy 


( 26 ) 
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As y —*■ oo it is required that 


v(y\x,t) 

dv(y,x,t) 


dy 


y—-oo 


y—i ►oo 


0 

0 


We use the generalized double Fourier transform defined as 


(за) 

(зб) 


/ o° r OO 

dt / dxv(y,x,t)e~ t( - (XX ~ wi (4) 

-oo J — oo 

where a is the complex wave number and u is the complex frequency. The lower limit 
of the time integration can be replace by 0, because v is a causal function in t ( v = 0 
for t < 0). The Fourier transform on time is then identical to the Laplace transform on 
time, with the transform variable S related to u> by the relation S = —iui. Thus (4) may 
also be referred to as a Lap lace- Fourier transform. 

Transforming Equation (1) yields the Orr-Sommerfeld equation 


{(D 2 - a 2 )(D 2 - A 2 ) + ietRU"} {$} = 0 (5) 

where 


A 2 = a 2 + iaR [t%) - ( 6 ) 



$ = 3>(y;a;,u;) 


$ is considered a function of y, with a and w as parameters. R is also a parameter, but 

omitted from the list of independent variables throughout this paper. 

Transformation of (2), the boundary conditions at the wall (y = 0), yields 

iu> 

$(0;a,w) = 


,2 — / j? 


a$(0;q,a>) _ Q 


dy 


(7a) 

(7b) 


7 



Transforming (3), the boundary conditions as y — ♦ oo , results in the requirements 



(8a) 

( 86 ) 


Equations (5) to (8) form a non-homogeneous system, where the nonhomogeneity appears 
in the boundary condition (7a). 

The solution of the equation is 

4 

* = ( 9 ) 

;'= i 

where a, to ) , (j = 1, . . . ,4), are the four fundamental solutions of the Orr- 

Sommerfeld equation. Cj = Cj(a,u) are constants with respect to y, to be determined 
by applying the boundary conditions. 

The outline of the formulation to this point is similar to G aster’s (1965), except that the 
latter is done in terms of the disturbance stream function (see Appendix C and Section 
6 ). 
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3. Solution in the transformed domain 

The fundamental solutions of the Orr-Sommerfeld equation need to be used in order 
to find the constants Cj from the boundary conditions. However, exact closed-form 
solutions of the equation for a general profile U (y) are not known, except for its asymptotic 
solutions as y — > oo. These solutions will be used when applying the boundary conditions 
(8) , and $ will then be expressed in terms of the fundamental solutions and their 
derivatives at the wall as follows. 


3.1. Asymptotic solutions of the Orr-Sommerfeld equation as y — ► oo 


As y — ► oo, U -* Uoo and u" 0. Taking U x = 1, the Orr-Sommerfeld equation (5) 
takes the form 

{(D 2 - a 2 )(D 2 - y 2 )}{$} = 0 (10) 

where 


y 2 = a 2 + iR(a — uj) 


( 11 ) 


and 



y — -oo 


2 - 2 


( 12 ) 

(123) 


Equation (10) is the asymptotic form as y — ► oo of the Orr-Sommerfeld equation, and 
can be solved exactly. Its four solutions are the asymptotic forms as y -+ oo of the 
four solutions of (5) , and are of the exponential form 


$i=e r * (j = l,...,4) 


(14) 


where rj are the four solutions of the characteristic equation 


(r 2 — ar 2 )(r 2 — y 2 ) = 0 


(15) 
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Let ri and r$ be the solutions of r 2 — fi 2 = 0 , and r-i and r 4 be the solutions 
of r 2 — or 2 = 0. Therefore solving for the rj 's requires taking the square root 
of the complex functions a 2 and /x 2 . To emphasize this point, we write r ; - as 


7*3,1 = ±(/X 2 ) 1/2 

(16) 

7* 4 ,2 = ±(or 2 ) 1/2 

( 17 ) 


Evaluation of the complex square root requires compliance with a convention for making 
the square root function single valued, and construction of branch-cuts in the a and ui 
planes. It seems that this point, although elementary, has been somewhat overlooked in 
prior work. The evaluation is detailed in Appendix A, where use of a mapping to an 
auxiliary plane yields clear definitions of the branch-cuts and a demonstration of their 
properties. It is shown there that 

7 * 3,1 = ±/x (18) 

where \i is the complex square root of (11), properly defined as a single valued function 
by the convention specified in Appendix A. 

The function /x depends on the three parametes a , w, and R. As shown in Appendix A, 
at fixed u and R , the branch-cuts for /z in the a-plane take the form of two hyperbola 

sections (Figures 2b and 6). In the w-plane , at fixed a and R , the branch-cut for /z is 

a straight line section (Figures 2a and 7). The choice of the branch-cuts in either case 
ensures that 

Real ( fj , ) > 0 for all a and u (19) 

T his property will be exploited when applying the boundary conditions. 

To assist in the evaluation of r 4) 2 , we define the function C as 

< = [(c* - ie)(a -He)] 1 ' 2 (20) 

where the convention for the complex square root function is in effect. It is shown in 


10 



Appendix A that 


( 21 ) 


Let 


r 4;2 = ± limC 
£—>•0 


a = a r + ioti 


( 22 ) 


where the subscript r indicates the real part of a, and i indicates its imaginary part. 
This notation applies also to other complex variables throughout this paper. Then (21) 
can also be written as 


r 4;2 = ±sign(a r ) • a 


(23) 


The branch-cuts which make C single valued are sections of the imaginary a- axis, from 
the branch-points ±ie to ±ioo (Figures 2b and 8). As e — *• 0 the branch-points approach 
the origin, and the branch-cuts span the whole imaginary a-axis, excluding the point of 
origin a = 0 , which stays a regular point. By the choice of these branch-cuts it is ensured 

that 

Real(C) > 0 for all a (24) 


This property will also be exploited when applying the boundary conditions. 

The asymptotic solutions are obtained by inserting rj from (18) and (21) into (14), 


= e" w 

(25a) 

$ 3 = e +w 

(256) 

$2 = lim e -( * y 

(25c) 

£ — 0 


$ 4 = lim e + ^ y 
£ — 0 

(25d) 


3.2. Application of the Boundary Conditions 

To apply the boundary conditions at y — * oo let p. = /x r + ifii and write (25b) as 

$ 3 = e^ y e^ iy ' (26) 
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The second exponential in (26) represents a pure oscillation, and the first is its amplitude 
modifier. Since y > 0 , growth or decay are determined by the sign of /i r . Since 
fi r > 0 for all a and u>, as follows from (19), the solution $ 3 is unbounded as y -> oo, 
requiring C 3 = 0 . By using a similar argument, following (24) , it is concluded that $i 
is unbounded as y — ► oo, requiring C 4 = 0 . $i and <J> 2 decay as y -*■ oo, therefore <J>i 
and $2 are acceptable solutions. 

The two remaining constants C\ and C 2 are determined by applying the boundary 
conditions (8) at the wall. The result for $ is 


where 


iu 

A 

Ao 

(27) 

A & A (y; a, w) = 4 >i$ 2 0 “ ^io ^2 


(28a) 

Ao = A(0;a,w) = $io$ 2 o - ^ 10^20 

(A at y = 0) 

(286) 

(j = 1,2) 


(28c) 

^•o = ^(0; a, w) (i = 1,2) 

($j at y = 0) 

( 28d) 

III 


(28e) 


It is important to emphasize the role of the particular choice of branch-cuts in this process. 
If different branch-cuts would have been constructed, then the growth and decay of the 
solutions as y — ► oo would depend on subdomains in the a and u planes, making the 
solution more complicated, if not impossible. 


12 



4. Inversion to the physical domain 


f.l. The singularities of $ 

The function $ posesses continuous and discrete singularities which are important for 
the inversion to the physical domain. The continuous singularities are the branch-cuts in 
the complex a and u planes which are discussed in in Section 3.1 above and in Appendix 
A. The discrete singularities are the zeroes of the two terms in the denominator of (27): 
(i) The zeroes of (w 2 - o»o 2 ) yield two poles at w = ±w 0 along the real axis of the the 
u>-plane; and (ii) the zeroes of A 0 yield additional poles in both the w-plane and in the 
a-plane. These latter pole singularities (ii) are discussed in the remainder of this section. 
These discrete singularities (ii) are obtained by solving 

Ao(a,w)=0 (29) 

For fixed R the expression (29) yields 

a = ctj (co) ( j = 1 ,■■■ ,N) (30) 

The expression (30) represents N w-dependent poles in the a-plane . As u varies, the 
poles trace trajectories there. Alternatively, (29), can be solved for u to yield 

io=uj(a) (jf = l,...,M) (31) 

which corresponds to M a-dependent poles in the w-plane, which form trajectories in the 
u>-plane as a varies. 

Both a and u> are complex in the expression (29), which can be considered as the 
generalized form of the dispersion relations which are obtained in the spatial or temporal 
normal mode analysis of linear stability theory. The eigenvalues in the spatial formulation 
are a special case of (30) with real u, and the eigenvalues in the temporal formulation 
are a special case of (31) with real a. 
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A complete description of the discrete singularities is needed in order to perform 
the inversion. The suitable form of this description are fixed Reynolds-number maps of 
families of constant u, and u> r trajectories of Oj(a;) in the a-plane, for all the j modes; 
or, alternatively, similar maps of constant a, and a r trajectories of ujj (a) in the w-plane. 

It is to be noted that negative frequencies and wave numbers need also to be included. 
It is, however, sufficient to specify in the a-plane only trajectories of a/(w) for w ’s with 
positive real frequencies u r > 0, because it follows from the symmetry property (B4) 
(Appendix B) that constant w,- trajectories of for u r > 0 are symmetric, with 

respect to the imaginary a-axis, to the ones for u r < 0. Similarly, it follows from (B5) 
that constant a, trajectories of Wj( a) for a r > 0 are symmetric, with respect to the 
imaginary w-axis, to the ones for a r < 0, and it is therefore sufficient to specify in the 
w-plane only trajectories of u>j(a) for positive wave numbers. 

We attempted constructing the required trajectories for boundary layer flow based on 
data from the open literature, but found that the available data are not sufficient for 
that purpose. The only reported computations known to us when both or and u are 
complex are by Gaster & Jordinson (1975) and Koch (1986). Some trajectories atj{u) 
can be derived from the data presented in those works, but unfortunately they are limited 
to negative values of u>i, a range which is not useful for the present analysis. The rest of 
the reported computations for boundary layers are of either the spatial or the temporal 
modes, where either or or a; is pure real. These computations (e. g. Mack (1976) ) show 
that for finite R the number of eigenmodes N and M is finite, and that there is only 
one unstable mode. The only theoretical proof is by Miklavcic & Williams (1982) and 
Miklavcic (1983) for the temporal case at finite R. This problem is not addressed here 
and we assume that N and M are finite. (1982) and Miklavcic (1983) for the temporal 
case at finite R. The problem of the number of modes is not addressed here and we 
assume that N and M are finite. For clarity we designate this mode with subscript j = 1 
and refer to it as the unstable Tollmien-Schlichting (TS) mode. We attach the subscripts 
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j = 2, . . . ,N (or M ) to the other modes, and refer to them as the higher, or stable, 
modes. 

Although the TS mode itself has been extensively investigated computationally, there is 
still lack of sufficient data to obtain the complete trajectory maps even for the spatial 
and temporal cases, because much less detailed data exist for the higher modes. 

It is, however, known that the trajectories of the higher modes u>j(a r ) and oij (ui r ) 
(j = 2,...,N or M ) , are in the lower half w-plane and in the upper half a-plane, 
respectively, for all Ft. In comparison, accurate trajectories of fti(wr ) (^* ^)> of 

t*>i(<*r) (<*«' = 0)> can be constructed for a wide range of w r and a r , respectively. At 
R > R C r the trajectory of wi(a r ) crosses the real w-axis and has a maximum in the upper 
half w-plane, and the trajectory of ai(w r ) crosses the real a-axis and has a minimum in 
the lower half a-plane. 

Computation of the generalized trajectories, when both a and ui are complex is 
beyond the scope of the present work, however, a qualitative investigation of the relative 
location of these trajectories for the unstable TS mode was performed in Ashpis h 
Reshotko (1986), based on the known spatial and temporal cases. It is shown there that 
the constant w t trajectory of ai(w) shifts upward with increased w,-, and is completely 
above the real a-axis for some positive Wj = cr . Similar qualitative conclusions regarding 
the trajectories of the higher modes cannot be drawn due to the unavailability of extensive 
computational results . In particular it is not known if trajectories can coalesce. 

Attention needs to be given also to the left-half complex planes, where the poles have 
negative real parts. The existence of eigenvalues with negative real parts, corresponding 
to waves moving in the upstream direction, has never been reported in literature. These 
are spatial eigenvalues a,-(w r ), fo: r positive u r , corresponding to poles in the left half 
a-plane, and temporal eigenvalues Wj(a r ), for positive a r , corresponding to poles in 
the left half w-plane. It is not clear if researchers have searched the left half-planes for 
eigenvalues without success, or if this search was not attempted at all. The implication 
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to the present work is that the existence of trajectories for u> r > 0 in the left half a-plane, 
and the existence of trajectories for a r > 0 in the left half w-plane is unknown. 
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4-2. Inversion formula 


Inversion of $ to the physical domain is obtained according to the inversion formula 

«(»;*,<) = tJ-v f f ^ < 32 ) 

JlJf 

where F is the inversion contour in the a-plane, and L is the inversion contour in the 
u>-plane. The inversion contours lie in the region of analyticity of $ in the respective 
planes, as shown in Figure 2, which incorporates results of the discussion of the previous 

section. 

Based on causality, the region of analyticity is an upper half-plane in the w-plane, with 
all singularities of $ lying in the remaining lower half-plane. Therefore L is a line passing 
above all the singularities of $ in the w-plane. 

Since v is defined for -oo > x > oo there is a strip of analyticity in the a-plane which 
includes the real a-axis. The strip is confined between ±ie , and all of the singularities 
of $ are located above and below it. As e — * 0 the strip degenerates to the real a-axis, 
allowed by the fact that the origin remains a regular point. Therefore F can be taken 
along the real a-axis, conveniently making a in (32) pure real. 

The two inversions in (32) are interelated. It means that if the w inversion is made 
first, a serves as a parameter whose value is taken along F. Since F coincides with the 
real a-axis, the singularities of interest in the w- plane are the ones for pure real a. The 
discrete ones are by definition the temporal eigenvalues, and the continuous ones are 
along the u>’s defined by the straight-line branch-cut for a = a r , as shown in Figure 2a. 
Since this branch-cut is always below the real u>-axis (Appendix A), L lies in the upper 
half w-plane above the maximum point of the temporal TS trajectory w x (a r ) (a,- = 0). If 
L is the line u> = ic , then a is positive. This conforms to the analyticity properties in the 
w-plane which were discussed above. If the a inversion is made first, the interdependence 
of the two inversions requires the singularities in the a-plane to be taken for values of u> 
on L, meaning for u>’s with = < 7 . The previous discussion showed that for these w’s 
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the poles are located in the upper half a-plane, and that the branch-cuts do not cross 
the real a- axis (Appendix A), which conforms to the requirement for existence of a strip 
of convergence there. 

At this point the inversion contours are determined and the inversion according to (32) 
can, at least in principle, be performed. 

J.3. The time- asymptotic form 

The present interest is principally in the time-asymtotic form of the solution, v as t — * oo. 
For this purpose we apply the method of Briggs (Briggs 1964) (see Section 1). When 
using this method, the inversion is made first from a to x, to obtain v as 

v(y;x,u) = J^(y; a, u)e t0ix da (33) 

followed by inverting v from u to t according to 



The method of Briggs requires deflecting the contour L towards L\, which is located 
slightly below the real w-axis (Figure 3a). The exponential term in (34) will make the 
integrand vanish along the straight portions of L\ as t — ► oo, leaving the time-asymptotic 
form to be determined by the singularities of v(y] x,u>) above L\. Simultaneously the 
contour F in the a-plane has to be deflected around constant u r trajectories of aj(u>) 
which cross the real a-axis as Uj — > 0. The trajectory of the TS mode displays this 
behaviour, therefore F is deflected to F\ , as shown in Figure 3b. The branch-cuts in 
the a-plane do not interfere with this process, because as was shown in Appendix A, 
the upper hyperbola branch-cut does not cross the real a-axis as long as w, > 0, the 
lower hyperbola branch-cut is below the real a-axis for all uj, and the imaginary axis 
branch-cuts do not depend on u at all. 
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Finally, the Briggs method requires one to check the possibility of the coalescing of 
constant w r trajectories originating at opposite sides of the real a-axis. Such an 
occurrence would be sufficient for the existence of a branch-cut of v in the w-plane, 
which corresponds to an absolute instability. Investigation of this possibility requires 
knowledge of all the trajectories a ; (w). As discussed in Section 4.1 above, from spatial 
stability calculations it is only known that for w, = 0 all of the higher modes are above 
the real a-axis. The possibility that their trajectories move into the lower half plane as 
Ui is decreased from its value on L, and coalesce among themselves or with the trajectory 
of the TS mode, cannot be excluded a priori, without performing extensive numerical 
calculations. This is left for future work, and at the present we make the assumption 
that coalescing of poles does not occur. The fact that absolute instabilities in boundary 
layer flow were never observed experimentally supports this assumption hueristically (it 
was not shown by the Briggs method that an absolute instability is sufficient for coalescing 

of trajectories). 

This assumption leads to the conclusion that v has no singularities above the real u>-axis, 
and L can then be deflected to L\ with the proper deflection of the contour around the 

poles ±u;o- 


Inversion from a to x 

First the inversion from a to x is performed according to 

v(y]x,u) = ^- [ $(y,cx,u})e iax da > (35) 

2tt J Fl 

To evaluate the integral we construct closed contours in the a-plane and apply the residue 
theorem. Two different contours are used, one for the domain x > 0, downstream of the 
ribbon, and the other for x < 0, upstream, as shown in Figure 3b. 

For x > 0 the contour is closed in the upper half a-plane with the semicircle Ti of radius 
r, deflected around the two branch-cuts. The closed contour consists of the section of 
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F\ from — r to -fr, of the sections of the semicircle Ti, and of the four sides of the 
branch-cuts, marked with the Roman letters I to IV (Figure 3b). 

The residue theorem yields 


/ + /+ /+/+/ +/ = 2iri Y. fi»Wy;“. w )« ,a l =et w (w) 

J-r J r x Jl JII Jin JlV j = 1 J v 


(36) 


where the integrand has been omitted, and AK u ) is the number of poles above F\. The 
first integral represents integration along the section of Fi from — r to +r. The rest of 
the integrals represent integration on the sides of the branch-cuts, and the the right hand 
side is the sum over the residues at the N discrete poles aj(u>), designated with the 
superscript (u). 

For the domain x < 0 a similar semicircle is constructed in the lower half o-plane, 
deflected around the branch-cuts there. Similar application of the residue theorem yields 

+r 

/ +/+/ + /+ / +/ =-2ii^R«s[%a,w)e ,ffI ] (i=tt ( ( ) w (37) 

J-r J r 2 Jv Jvi Jvil JviII j =1 1 

The roman letters designate the sides of the branch-cuts as marked in Figure 3b and the 

residues are calculated at the poles below Fi , designated with the superscript (£). 

As r —* oo the first integral in the the left hand side of (36) and (37) approaches the 
integral in (35) . The integral on Ti vanishes for x > 0 as r — * oo, and the integral on T 2 
vanishes for x < 0 as r — » oo. Assuming all poles are of first order the following result is 
obtained 

v(y,x,uj ) = v D (y,x,u)) + v c {y,x,v) (38) 

wherein vq is the discrete spectrum in the u-plane, given as 


VD(y\x,u) = 

jv(«> 

A(y;aj u) (w)>“) 

x > 0 

(39a) 

0 

dAo(a$ w) (w),w) 





da 



VD(y,x,u) = 

N w 

y > U) 

A(y,a^\u),w) i<K (f) (u)m 

x < 0 

(396) 

3 = 1 0 

dA 0 (a^\u),uj) 


da 
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and vc are the continuous spectra 


v c (y,x,aj) = vc 1 (y-,x,u >) + vc?(y,x,ui) (40) 

where subscript 1 designates spectra originating from integration around the imaginary 
axis branch-cut, and 2 designates the spectra originating from the hyperbola branch-cuts 
as follows: 

x : u) = — — J $(y; a,u>)e tax da — — $(y; a,u)e da, 

x > 0 (41a) 

0c 2 (y; x,u>) = — — J $(y,a,u>)e ,ax da - — j^^{y,Oi,oj)e iax da , 

x > 0 (416) 

6 Cl (y;x,u>) = J $(y; a, u)e tax da - ^ $(y, a, u)e' ax da , 

x < 0 (42a) 

**(»:«■") = ~i ' 

X < 0 (426) 


If. 5. Inversion from ui to t 

Next v is inverted from w to < according to 

v(y\x,t)= f v(y,x,uj)e~ lut cL > ( 43 ) 

JLi 

The integral vanishes on the straight section of L\ for t — ► oo, leaving v to be calculated 
by integrating around the circles surrounding the poles ±wo- Using the residue theorem, 

v(y-x,t) = 5^Res[w(y;x,w)e _ ‘ a ’ t ] u)=;two ( 44 ) 
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Detailed evaluation of (44) is given in Ashpis k Reshotko (1986) and is outlined in the 
following. 

By observing (39) , (41) and (42) , it is seen that v consists of sums and integrals of 
functions of the form 

G = 2 g(y; x, w) (45) 

u 2 -UJ$ 

It is found that 

Res [G] w=±u> 0 = ± ^(y; ±w o) (46) 

Using the symmetry property of g, 

g(y ; x,u 0 ) = g(y, x, -w 0 ) (47) 

where the overbar represents the complex conjugate, it is obtained that 

Res [<?] w=±Uq = Real[g(y, x, w 0 )] (48) 

Applying these results in (44) leads to the final result for the time-asymptotic form of v 
as follows 


v(y\x,t) = v D (y\x,t) + v c (y,x,t) (49a) 

vc(y,x,t) = v Cl {y,x,t) + vc 2 (y,x,t) (496) 


where the various terms are detailed in the following. The added superscripts (u) and 
(£) designate the downstream (x > 0) and the upstream (x < 0) regions, respectively. 
The discrete part of v, designated vp, is: 


VD\y\x,t) 


N^) 

-E Im 

i= i 


i 


A(y; q( u) (u;o),u>o) ^ ra («) (a , 0 ' )r . h , ofl 
6A 0 (a( u) (u;o),wo) 

da 


> , x > 0 (50a) 




nW 


A(y;af\uo),u 0 ) Jfa (< ) (wo)x _ Uot] 
dAo(af\u>o),uo) 


j = i 


e 1 j 


da 


x < 0 (506) 


22 



where the sign designates that this is the time-asymtotic limit. A non-zero value for 
(50b) is contingent on existence of poles in the left half-plane as discussed before. 

The continuous part of v originating from the imaginary-axis branch-cut, designated vc x , 


is: 


~ Im 
27r Im 



A ( i l \y,<r,u;o)e cx dcr 
A ( / ) (y;<7,u; 0 )e +<rr da 



where 

.( u )/ x A(y; i<T+^j 0 ) A(y; i<r~, cj 0 ) 

i (y.^.wo) a 0 ( i<r+,u o) Ao( i<r - ,u;o) 

.(0/ x A(y; — icr + ,u) 0 ) A(y;-»<r~,uJo) 

i (y.ff.wo) Ao(— i<T + ,wo) Ao(— i<r - ,uJo) 

wherein 


x > 0 (51a) 
x < 0 (516) 

(52a) 

(526) 


zcr* = lim (e + icr) 

£^0± V 

The ± signs correspond to the sides of the imaginary axis branch-cuts, as shown in Figure 


3b. 

The continuous part of v originating from the hyperbola branch-cuts, designated vc 2 , is 


v$(y,x,t) ~ A Real /,-•>(«)* jr4“ ) (»i<r. U0 )«-”e i[, * (w) *-"* ,1 <to} , 

x > 0 (53a) 

„<§(»;,,*) ~ A Real jT d<r} , 

X < 0 (536) 


where 


x ra«i(<r,w 0 ) 
^2 ; (y;^.wo) = y — Yg — 


+ i 


A(y; <*£,(*, m)),m>) 

A o(aB 1 Kw 0 ),uo) 


A(y;Q Bl (a-,^o),^o) 

Ao(otB 1 (o’» w o),w 0 ) 


(54a) 
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A%\y;<r,u 0 ) = 


862(0, w q ) 

da 


A (y; <XB 3 ( tr > u o)> u o) _ A (y? 

Ao(ar£ 2 ( 0 -,w o ),w o ) Ao(4 2 (<r, w 0 ), wo) 


( 546 ) 


cr is an integration variable, subscript B\ designates the upper hyperbola branch-cut, 
Bi the lower one, the superscripts ± designate the side of the respective branch-cut, as 
marked in Figure 3 b. 

612 and ai j2 are the coordinates of the branch-points ctBP lt3 — (^1,2 + ^1,2) for w = wo, 
and are given by the following equations 

1 / 2 ' 


R 

2(^0) = y { t 1 + 


H 1+ H 


4w 0 

~R 


(55) 


6 l, 2(^0) = 


W 0 


1 ±2 


O~l, 2 (^o) 

R 


W 0 


For — << 1 it is obtained that 

R 


(Tl(wo) « ^ 
u? 

61 (wo) ~ w 0 
62 (wo) ~ -Wo 


(56) 


(57a) 
(576) 
(57c) 
(57 d) 
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5. Description of the results 


The expressions (50) to (56) show that the vibrating ribbon excites discrete and 
continuous spectra, whose frequency in the time-asymptotic limit is equal to the frequency 
of the ribbon. 


The discrete spectrum consists of the spatial eigenmodes of the flow. These are 
discrete travelling waves with frequency wo and wavenumber Qfj r (wo). The phase velocity 
of the wave is 


The growth factor is e a ii( u °) x 
which the mode is excited, is 


Cj = 


Wo 


(58) 


a,, (wo) 

and the coupling coefficient, indicating the extent to 


_ A(y;o j (wo),w 0 ) (59) 

J dA(aj(wo),w 0 ) 

da 

The quadrant in which the pole is located then determines the properties of its 
corresponding wave. Its region of influence is determined by its position above or below 
the contour F x . The various possibilities are depicted in Figure 4. Poles in the right 
half q-plane correspond to waves propagating in the positive x direction (positive phase 
velocity), and the ones in the left half q-plane correspond to waves travelling in the 
negative x direction. Growth in the region x > 0, downstream of the ribbon, is if the 
corresponding pole is below the real q-axis; and in the upstream region , x < 0, growth 
in the negative x direction is if the pole is above the real q-axis. The right half q-plane 
was thoroughly investigated numerically by various researchers (see discussion in Section 
4 above), and the cases described in Figure 4 (a-ii) and (d-i) where never found. The 
cases (b-ii) and (c-ii) cannot be expected to be found when extending these computations 
to the left half q-plane. 

For x — » oo, the higher modes decay spatially, leaving the Tollmien-Schlichting mode 
to be the dominant mode in the downstream region, and the summation in (50a) can be 
omitted. 
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The continuous spectra affect the regions upstream and downstream of the ribbon. In 
each region it consists of two parts; the one described by (51), consisting of a superposition 
of standing waves, and the one given by (53) is a superposition of travelling waves. In 
both cases the waves decay away from the ribbon asi-> ±oo, although growth is possible 
for finite x close to the origin. 
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6. Discussion 


In the work presented here we have revisited the vibrating ribbon problem which 
was first done by Gaster (1965). For the purpose of comparing the two works, we use 
the notation G when referring to G aster’s work, and present its outline in Appendix C. 
We have written the equations and the boundary conditions in terms of the disturbance 
normal velocity, while G used the disturbance stream function. As shown in Appendix 
C, the two formulations are equivalent and lead to the same equation and boundary 
conditions. Except for this minor difference, we followed the formulation of G up 
to the point where the double transform is applied (Equation (9) ). From there on 
our solution method differs considerably. We expressed the solution in terms of the 
fundamental solutions of the Orr-Sommerfeld equation, exercised some care in the 
evaluation procedure, and applied the Briggs method to obtain the time-asymptotic 
solution. In G the solution is written in terms of a nonspecified complete solution of 
the the Orr-Sommerfeld equation and boundary conditions , and the inversion contour is 
not determined in a straightforward way (see Appendix C for details). 

A significant outcome achieved by the present work is that the solution includes the 
continuous spectra, which are missing from G. The wave components of these spectra 
are identical to the spatial continuous spectra obtained in Grosch and Salwen (1978) 
and Salwen and Grosch (1981) by looking for pure oscillatory eigenmodes of the the 
Orr-Sommerfeld equation . Here these eigenmodes emerge in a natural way, simply by 
constructing branch-cuts in the complex planes, as required whenever a complex square- 
root function is encountered. 

The discrete spectrum obtained in the present work is equivalent to the one in G, 
with the exception of the factor of 1/2 which appears there ( compare (Cll) to (50a) ). In 
addition, the result in G depends on assumptions involving the sign of the group velocity 
(see Appendix C). No such asssumptions are used in the present analysis. The reason 
for the discrepancy by the factor of 1/2 is an error in the computation of the residue at 
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u> — u)q: The inversion contour in G was taken along the real w— axis (Figure 9). In that 
case the following expression, appearing in (C6) takes the form 


lim — 2 = 2 r 2 - *x[^(Wr + w 0 ) + 6(io r - w 0 )] (60) 

w, — 0 U) z — UIq Wj? — Uq 2 

The contribution from the term containing the 8 functions, which is equal to the 
contribution from term preceeding it, was not taken into account in G. 

It was obtained here that the ribbon has a downstream influence by excitation of 
both the discrete and continuous spectra. The new finding is that it also has an upstream 
influence via the continuous spectra . It was shown that the possibility of an additional 
upstream influence by excitation of the discrete spectra needs also to be considered. This 
can be verified or excluded only by numerical search for poles in the left half plane, which 
correspond to spatial eigenvalues with negative real parts for positive frequencies. The 
spatial decay in the time-asymptotic limit of the continuous spectra is much greater in 
the upstream direction compared to the downstream direction. This is concluded from 
(57), showing that at large Reynolds numbers a 2 » ci, which affects the exponential 
decay terms preceding the integrals in (53). 

The distinction between singularities affecting the upstream region and the ones 
affecting the downstream region is useful also in clarifying the difficulty expressed in 
Salwen and Grosch (1981). They point out that the wave component of the continuous 
spectra which propagates upstream from x — 00 and the standing waves whose amplitudes 
increase towards x = 00 are not physically acceptable. This difficulty is resolved by the 
analysis here, which shows that these waves are excited only in the upstream region of 
the ribbon and therefore are physically acceptable. 

Once solved systematically and correctly, the complex plane contains all the discrete 
and continuous singularities, and the the solution can serve as proof of completeness of 
the eigenfunction expansion in the spatial case. This is in analogy with the proof of 
completeness for the temporal case done by Salwen Grosch (1981) by comparing to the 
solution of Gustavsson (1979). Although proof of completeness for this case is reported in 
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an abstract by Salwen, Kelly and Grosch (1980), details were never published and could 
not be obtained from the authors. 

Our solution also sheds light on the difficulties in the IBVP solutions reviewed in 
Section 1. In these solutions a Laplace transform was applied to the longitudinal spatial 
coordinate x, which is equivalent to the Fourier transform of a function causal in x (A 90 
degree rotation relates the Laplace 5-plane to the Fourier a-plane). Therefore, formally, 
the application of the Laplace transform should yield the same discrete and continuous 
singularities as the Fourier transform. However, the complex 5-planes in the work of 
Tsuge k Rogler (1983) and of Aldoss (1982) do not include the counterparts of the 
imaginary axis branch-cuts. In addition there is an inconsistency in these works, as well 
as in Tumin k Fedorov (1984), where branch-cut singularities protrude into what should 
be a half plane of analyticity of the 5-plane. As is shown, the disturbance introduced by 
the vibrating ribbon has an upstream influence, at least through the continuous spectra, 
and so will any disturbance introduced at x = 0. By applying the Laplace transform this 
possibility is excluded, and it is suggested that this is the source of the mathematical 
inconsistency, and that it is incorrect to use the Laplace transform in these problems. 
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Appendix A. The functions fi and ( and evaluation of rj 

The essentials of the evaluation of the functions // and £ and the four solutions rj of the 
characteristic equation (15) are contained in this Appendix . The detailed evaluation can 
be found in Ashpis & Reshotko (1986). 

A.l The complex square root function 

rj are evaluated as the complex square roots given in Section 3.1 as 

r-3,i = ±(^ 2 ) 1/2 (16) 

r 4j 2 = ±(a 2 ) 1/2 (17) 

We follow in detail the well known elementary steps to evaluate the square roots. To 
excercise this extra care we first review the complex square root function F = Z 1 / 2 . F 
has two single valued branches F + and F ~ , which are related by F + = —F~. To make F 
single valued, it is needed (i) to construct a branch-cut in the Z-plane, (ii) to characterize 
F + and F~ by a specified definition, and (iii) to select one of the branches F + or F~ to 
represent F. 

These three steps are executed as follows: (i) Any line from the branch-point Z = 0 to 
oo is an acceptable branch-cut, and we select the negative real Z-axis as the branch-cut. 
(ii) We define F + to be the branch of F that yields real positive F at real positive Z, F~ 
is then the branch that yields real negative F at real positive Z. (iii) Finally, we choose 
the single valued function to be represented by the branch F + . This process establishes 
a convention for defining the single-valued complex square-root function which is used 
throughout this paper. 

This convention ensures that 

Real(F) > 0 for all Z (Al) 
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The situation is illustrated in Figure 5 . It is easy to demonstrate that (Al) will not be 
valid for any other selection of a branch-cut in the Z- plane. 


A. 2 evaluation of r\ and 7*3 


Let Z be defined as 


Z = a 2 + iR(a — u) 


(A2) 


then, from (11) 


p = Z l/ 2 


(A3) 


where p is a single-valued function, according to the convention described above. 
Following (16) it is obtained that 

7*1 = —p (A4) 

7*3 = P (^ 5 ) 


Z depends on or, u and R. For fixed R, either u, or a , can be taken as a parameter. 


A. 2.1 p with uj as parameter 

The case where u is the parameter is discussed first. In this case p can be viewed as a 
mapping from the cv-plane to the /i-plane via the auxiliary plane Z . The various planes 
are illustrated in Figure 6. 

The mapping from the a-plane to the Z-plane is investigated in the following. 

Let 


Z — Zr “1" iZi 

(AS) 

a = ot r + i(*i 

(Al) 

U) =■ U> r + iu>i 

(AS) 
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Substitution into (A2) and equating the real and imaginary parts yields 



Z± -|- Ridr 
2a r 

d2 

f±- + Rui-Zr 
4 


(A9a) 

(A9b) 


which shows that constant Z T and Z,- lines are mapped into hyperbolas in the a-plane. 
(A9a) represents a family of canonical hyperbolas of the first kind, and (A9b) represents 
a family of canonical hyperbolas of the second kind, both shifted to or, = —R/ 2. In 
addition to the independence, (A9a) depends only on u> r , the real part of u>. (A9b) 
depends only on its imaginary part u The mapping is illustrated in Figure 6 for the 
case of a pure real and positive u>. The real Z - axis is mapped to the hyperbolas a\ and 
a 2 ; the imaginary Z-axis is mapped to the hyperbolas 6i and 6 2 . The quadrants of the 
^-plane and the domains in the a-plane into which they are mapped are marked with the 
roman letters I to IV. The branch-point Z = 0 is mapped to the points BP\ and BPi- 
These points are the intersection points of the hyperbolas ai with b\, and of a 2 with f> 2 .It 
then follows that the branch-cut in the Z-plane is mapped into the two marked sections 
of the hyperbolas a\ and a 2 , from the branch-points BPij to ±oo. 

The coordinates of the branch-points and equations of the branch-cuts are given in Ashpis 
& Reshotko (1986). The branch-points vary with u> as follows: As w,- increases the 
branch-points move towards ±ioo along the hyperbolas a\ and a 2 , respectively. The 
upper branch-point BP\ is above the real a-axis for Wj > 0. As u>,- decreases to negative 
values BP\ moves below the real a- axis. The lower branch-point BP? stays below the 
real a-axis for all cj; . For u> r < 0 the picture is symmetric with respect to the imaginary 
a-axis. 

It follows from the previous discussion that the auxiliary Z-plane, with the branch-cut 
along the negative real axis is mapped into the right half //-plane. Since the whole the 
a-plane is mapped to the Z-plane, it follows that with the specified branch-cuts it is 
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ensured that 


Real (fi) > 0 for all a and w. 


(- 410 ) 


A practical formula for calculation of n is given by 


A* 




where 


Z r = a 2 - a 2 - R(ai - u>, ) 
Z { — ”1* (^r) 

r = (Z r 2 + z?y > 2 


(-411) 


(-412) 

(-413) 

(-414) 


The -f sign is used in regions I and II, the — sign in III and IV of the a-plane (Figure 6). 

A. 2.2 pt with a as parameter 

When a is taken as the parameter in (A2) with fixed R , // can be viewed as a mapping 
from the w-plane to the //-plane via the auxiliary plane Z , as depicted in Figure 7 . We 
obtain in a similar fashion 


u>i = — {Z r — o? r + af + Ron ) 
K 


(-415) 


= 4(-Zj + 2a iar + Ra r ) (-416) 

R 

These expressions show that constant Z r and Z x lines in the z- -plane are mapped into 
straight lines in the w-plane. 

The branch-point Z = 0 is mapped into the point BP 0 whose coordinates are 
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WBP 0 = *?o + i<r o 


(417) 


where 


Oo = -^[“r - (“? + R<*i)] 


m = a r 



C-A18) 

(- 419 ) 


The branch-cut extends from BPo to — oo in the lower half w-plane parallel to the 
imaginary axis. For a, = 0 the branch point forms the trajectory cr 0 = — ^7? 2 j therefore 
the branch-cut is always below the real w-axis. 

This branch-cut is mapped into the branch-cut in the Z-plane therefore it follows that 

Real (/i) > 0 for all u and a. (A20) 


A. 3 Evaluation of and the function £ 

A similar analysis is applied to (17) and (21). As noted in Section 3.1, it is useful to 
avoid ambiguity by defining the function £ as 


to 

II 

IT 

i 

IT 

+ 

(421) 

(For a similar idea see Carrier et al (1966) pg 315, 346). 


Then 


a 2 = lim £ 2 

f— 1-0 

(.422) 

We define z as 


z — {oi — ie)(a + is) 

(423) 

and it follows that 


II 

b* 

i-* 

to 

(424) 
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where £ is a single valued function subject to the previous convention regarding the 
complex square-root. Therefore 

r 2 = — lim C (A2b) 

£ — 0 

r 4 = lim C (A26) 

£ — 0 

Similar to the treatment of //, we view ( as a mapping of the a-plane, via the auxiliary 
plane z (Figure 8). The process here is, however, simpler because < does not depend on 

R or u>. 

As before, we select a branch-cut in the z-plane along the negative real axis, therefore 
the whole z-plane is mapped into the right half £-plane. 

The mapping from the a-plane to the z-plane is given by 


where 




(A27a) 

(A27b) 


z = z r + izi 


(A28) 


(A27a) shows that constant z r lines in the z-plane are mapped into a family of canonical 
hyperbolas of the second kind in the or-plane. (A27b) shows that constant z* lines are 
mapped into a family of canonical hyperbolas of the first kind . The branch-point z = 0 
is mapped to the points ±ie in the ct-plane, and the branch-cut is mapped to the sections 
of the imaginary a-axis, from the branch-points to ±oo, as shown in Figure 8. Selecting 
these sections as branch-cuts ensures that 


Real(C) > 0 for all a 


(A29) 


A working formula for the calculation of £ is given as 


c = 






a 



(A30) 
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where 


= a 2 r - of + e 2 

(A 31) 

= 2a r a,- 

(A32) 

= « + 

(A33) 


and the + is taken in the quadrants (1) and (3), and the — in qudrants (2) and (4) of 
the a-plane, as marked in Figure 8. 

As e — * 0 the two branch-points approach the origin a — 0. This point, however remains 
a regular point of (. Therefore in the limit e — ► 0 the branch-cuts coincide with the 
whole imaginary a- axis, with the exlusion of the point of origin. 

Also as e — *• 0 (20) can be written as 

£ = sign(a r ) ■ a (-^34) 
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Appendix B. Symmetry properties 


The function $ exhibits the following symmetry properties 

$(y;a,w) = 0(y;— a, — u7) (^ 1 ) 

where the overbar designates the complex conjugate. This property follows from the 
Orr-Sommerfeld equation and also from the properties of the double Fourier transform 

of a real function. 

The single Fourier transform of the real function v has the property 

v(y ; x,u) = v(y; x, -ZJ) ( B<2 ) 

It follows that 

Ao(ct, u) = A 0 (— -w) 

and that 

ttj(w) = -aj(-uJ) ( 54 ) 

Wj (a) = -wj(-a) ( 55 ) 

Additional details related to these properties are given in Ashpis & Reshotko (1986). 
It is demonstrated there that the symmetry properties are not valid with a choice of 
branch-cuts other then those specified in Appendix A. 
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Appendix C. Outline of Gaster’s solution (1965) 


Gaster’s formulation and solution of the vibrating ribbon problem ( Gaster (1965) ) is 
summarized in this appendix. For brevity we here refer to Gaster (1965) as G. 

We found some details in G somewhat difficult to follow, mainly due to to the lack 
of accompanying figures describing the complex planes. It was also noted in private 
communication (Gaster, 1984) that G contains some printing errors and has an omission 
of a short paragraph. For the purpose of comparison between G and the present work, we 
find it useful to outline the solution of G in this Appendix and to provide a supplementary 
figure. 

Let us designate the disturbance stream function as ip(y;x,t) and its double transform 
as \?(y; a, w), and use the notations of the present work: Then the disturbance stream 
function $ in G is ^ in our notation; the frequency /? in G is u in our notation; and the 
frequency of the ribbon u? in G is u>o in our notation. The other pertinent notations are 
identical. 

The formulation of the problem in the present work follows the one in G, except that in 
G the formulation is in terms of the disturbance stream function ip, while here it is done 
in terms of the disturbance normal velocity v. 

After a double transform, G obtains the Orr-Sommerfeld equation for 

(U(y) - w/c*)(tf" - a 2 *) - u"{y)V - (-i/aR){V IV - 2a 2 *" + a*9) (Cl) 
subject to the boundary conditions 
at y = 0: 

= <C2a) 

att(0;a,a;) _ Q (C26) 

oy 
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as y — > oo: 



T ^ n 

*• * 

(C 3a) 

since 

dV 

v ~ dx 

(C36) 

then 

$(y\a,io) = -ia^!(y,a,to) 

(C 4) 

From (C2a) 

to 1 

to 2 - uol ^(0; ct.w) 

(C5) 

which yields after substitution into (C4) 



ito tf(y; a,w) 
$(y,a,u) - w 2 _ w 2 *( 0 ;a )W ) 

(C6) 


The inversion to the physical domain in G differs from the one performed in the present 
work in the way the inversion contours are determined, and in the use of the symmetry 


of ^ with respect to a and to. 


Two possible inversion contours in the or-plane are considered by G, and the correct 


contour is selected based on results of the next inversion from to to t. The inversion from 
a to a; is essentially similar to the one done here, although G uses the symmetry property 
of with respect to real a and to, but the inversion from u to t is different than the one 
used here. The closed contour constructed by G for application of the residue theorem 
is shown in Figure 9. We drew this figure based on the description in G, and provide it 
here for illustration. It is argued in G that the inversion integral vanishes for t oo on 
the section A 2 of the contour. This argument utilizes assumptions which involve the sign 
of the group velocity. 

The result for the dominant TS mode obtained in G is 


v(y,x,t ) = -Im 


1 ^(y;a(a;o),a>o) c i[ 0r ^ 0 )-q, 0 tl > 

2 ^(Oja^o),^) 

da 


(C 7) 


39 



^ is a linear combination of the fundamental solutions of the Orr-Sommerfeld equation. 
The constants, corresponding to the Cj ’s of (9) in the present work, were not evaluated 
in G. 

In order to compare the results we write (C7) in terms of the fundamental solutions $j 
and A (see (9) and (28)), as follows: 

Using (C4) for the \P-ratio in (C6) yields 


^(0;a,w) $(0;a,u>) 

substitution of (27) into (C8) yields 

^(t/;a,^) _ A(y;or,u;) 
\k(0;a,w) A 0 (a,u;) 


then (C6) is 




iw A(y;a,u;) 
u 2 — u>q Ao(a,w) 


and (C7) takes the form 


(C 8) 


(C9) 


(CIO) 


v(y-x,t) 


I 1 A(y:a(^)^e) , 

lm j 2 ^Ap(a(uJo),a*o) e 

^ da J 


(Cll) 
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Figure 1. Model of the vibrating ribbon problem. 
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FIGURE 2. 


Figure 2. Undeflected contours for the inversion integrals (schematic). 

The points marked with o in (a) indicate poles for the <x r marked 
with □ in (b). The curves passing through these points are 
trajectories u>j(a r )• Points marked with x in (b) are poles for u 
marked with x in (a). The curves passing through these points 
are trajectories ocj{ut)- The branch cut in (a) is for a = a r . The 
branch cuts in (b) are for w = u>jj. The question-mark symbols ? 
in the left half-planes designate yet unknown existence of poles. 
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Figure 3. (a) Integration contour in the u;-plane. 

(b) Closed integration contours in the a-plane. 
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FIGURE U. 


Figure 4. Effects of the location of the pole in the a-plane on its 
corresponding wave for u> r > 0. 

x designates the ctj for u>j = 0, o designates ocj for to = to^ . 
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FIGURE 5. 


Figure 5. Illustration of the mapping F = for the branch F + . 

Roman letters indicate the quadrants in the Z-plane and their 

image in the F-plane. 
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(a) Cb) (c) 

FIGURE 6. 


Figure 6. The mapping from the a-plane to the ^-plane via the auxiliary 
plane Z. 

The hyperbolas a\ and a ,2 are mapped into Z{ = 0, the hyperbolas 
b\ and b 2 are mapped into Z T = 0. the (+) and (— ) indicate the 
two sides of the branch cut. Roman letters indicate the quadrants 
of the Z-plane and their images in the /z and a planes. 
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(3) (b) (c) 


FIGURE 7 . 

Figure 7. The mapping from the u?-plane to the /i- plane via the auxiliary 
Z- plane. 

The line a is mapped into Zi = 0, and the line b is mapped into 

Zr = 0. 
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(b) 

FIGURE 8. 


(c) 


Figure 8. The mapping from the a-plane to the £ -plane via the au xili ary 
plane z. 

The real and imaginary a axes are mapped into z,- = 0. The 
hyperbolas b\ and &2 are mapped into z T = 0. 
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FIGURE 9. 


Figure 9. Illustration of G aster’s (1965) inversion contour in the w-plane. 
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